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Cooperative Emission of a Pulse Train in an Optically Thick Scattering Medium 
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An optically thick cold atomic cloud emits a coherent flash of light in the forward direction when 
the phase of an incident probe field is abruptly changed. Because of cooperativity, the duration 
of this phenomena can be much shorter than the excited lifetime of a single atom. Repeating 
periodically the abrupt phase jump, we generate a train of pulses with short repetition time, high 
intensity contrast and high efficiency. In this regime, the emission is fully governed by cooperativity 
even if the cloud is dilute. 

PACS numbers: 42.50.Md, 42.25.Dd 


The seminal work on superradiance of R. Dicke in 1954 
has opened up tremendous interest in studying coopera¬ 
tive emission of electromagnetic radiation from an ensem¬ 
ble of radiative dipoles (see [T] for the original proposal, 
mig for reviews and mini for recent related works). In 
his original proposal, R. Dicke considered an ensemble 
of N excited two-level atoms confined inside a volume 
smaller than A^, where A is the transition wavelength. 
In this context, a macroscopic polarization is built up 
in the medium upon incoherent spontaneous emission. 
This Dicke superradiance mechanism leads to the co¬ 
herent emission of an intense pulse with a decay time, 
Tp) = that is shortened by a factor of N~^ 

with respect to the atomic excited state lifetime, T”^. 
For practical implementation in the optical domain, the 
Dicke model was extended to media with volume larger 
than A^ mnn [12]. In those cases, the propagation of 
the electromagnetic field in the medium and the spatial 
mode density must be taken into account. If the medium 
is dense, i.e., pX^ I, where p is the radiator spatial 
density, it still exhibits the main feature of the Dicke su¬ 
perradiance, namely, the emission of a short pulse after 
some delay m H im E]. It was, however, pointed out 
in El, that the superradiant pulse decay time should be 
corrected as t = Tp)/p. /r < I is a geometrical factor 
corresponding to the solid angle subtended by the super¬ 
radiant emission Eiiia. 

For a dilute scattering medium, i.e., pX^ <C I, the Dicke 
superradiance mechanism does not occur |I5) . Never¬ 
theless, an optically thick medium driven by a coherent 
incident field shares interesting similarities with Dicke 
superradiance; here, the cooperativity factor Np is re¬ 
placed by the optical thickness of the medium [niiisi. 
Once a driving coherent field is abruptly switched off, 
like in a free induction decay (FID) experiment [2(11127j . 
a short coherent cooperative flash of light is emitted in 
the forward direction. The flash duration is inversely pro¬ 


portional to the optical thickness and the bare linewidth 
of the transition |26j . A similar phenomenon occurs for 
the optical precursor, i.e., when the driving coherent field 
is abruptly switched on [55] ■ 

In a coherently driven medium, the incident probe fre¬ 
quency can be detuned with respect to the atomic reso¬ 
nance, leading to a nontrivial phase rotation of the coop¬ 
eratively emitted field (see [27] in the optical domain and 
[55H5T] for 7 -ray pulses in Mossbauer spectroscopy exper¬ 
iments) . In this Letter, we report the generation of high 
repetition rate and high intensity contrast pulse trains in 
an optically thick cold dilute atomic ensemble using the 
setup schematically shown in Fig.j^a). An example of a 
pulse train, generated in our experiment by periodically 
changing the probe phase, is shown in Fig. [^b). As a 
consequence of cooperative emission, the repetition time 
Tr of the pulse train can be shorter than the atomic ex¬ 
cited state lifetime, F”^. Moreover, we show that at high 
repetition rate, the single atom fluorescence is quenched. 
This constitutes a rather counterintuitive result where 
the emission in free space is fully governed by coopera¬ 
tivity, in contrast with the usual situations where it is 
enhanced by a cavity surrounding the medium [32]. 

The scattering medium is a cloud of laser-cooled ®®Sr 
atoms (see [33] and [27] for the details of the cold 
atoms production line). The ellipsoidal shape of the cold 
cloud has an axial radius of 240(10) p.m and an equa¬ 
torial radius of 380(30) qm, with peak density around 
4.6 X 10^^ cm“^ for a total of 2.5(5) x 10® atoms. A = 
689 nm is the wavelength associated to the t ^Pi in¬ 
tercombination line (bare linewidth of r/27r = 7.5 kHz) 
used in this experiment. pA® = 0.15, which puts us in 
the dilute regime. The temperature of the cold gas is 
T = 3.3(2) pK. We get kv = 3.4F, indicating a signifi¬ 
cant Doppler broadening of the narrow intercombination 
line, k = 2 tt/X is the wave vector of the transition, and 
V is the rms velocity of the gas. The optical thickness 
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Figure 1. (color online), (a) Experimental setup: a laser 
beam is sent through an acousto-optic modulator (AOM), 
which may be used to switch on and off the incident beam, fol¬ 
lowed by an electro-optic modulator (EOM), which abruptly 
changes the phase of the incident probe field, Eq. (b) A pulse 
train generated at a repetition time of Tn = 0.12r“^ by a 
periodic abrupt phase change of tt. Here, the probe laser 
is at resonance, (c) Electric fields just before and just af¬ 
ter an abrupt phase jump of n are represented schematically 
in the complex plane. Before the phase jump, the forward 
scattered field Es destructively interferes with Eq. After the 
phase jump, they constructively interfere. The transmitted 
field is denoted by Et. 


depends strongly on the temperature. We measure 19(2) 
along the equatorial axis at resonance. 

A 150 /im diameter probe laser beam, tuned around 
the intercombination line, is sent through the cold 
atomic gas along an equatorial axis. The probe 
power is 400(40) pW, corresponding to 0.45(5)/sat 
{Isat = 3 p,W/cm^). We measure the forward transmit¬ 
ted intensity of the probe using a photodetector, in¬ 
tegrating over the transverse dimensions of the trans¬ 
mitted beam. We apply a bias 1.4 G magnetic field 
along the beam polarization during the probing phase, 
making the atom an effective two-level system on the 
TO = 0 —f ^Pi,m = 0 transition. 

The ellipsoidal shape of the cloud is modeled by a slab 
geometry, so that the coherent transmitted electric field, 
in the frequency domain, is given by 


In the above equation, n(oj), Eq, c, and L are the complex 
effective refractive index, the incident optical field, the 
speed of light in vacuum, and the slab thickness along the 
laser beam, respectively. For a dilute medium, n(uj) = 
l + pa(Lu)/2 [33], with the two-level atomic polarizability. 


Q:(a;) 


SttFc^ 1 exp (-uV2f2) 

S-kv + tr/ 2 - 


6 = uj — ujo is the detuning of the probe laser frequency 
w with respect to the bare atomic resonance frequency, 
ujq. The effect of Doppler broadening is included in 
the polarizability by averaging over the thermal Gaus¬ 
sian distribution of the atomic velocity v along the beam 
propagation direction. The transmitted intensity It(t) is 
computed following [27j . and by performing an inverse 
Fourier transform. We define, for given 6 and v, the op¬ 
tical thickness by(S) and the relative phase 0g(i5) between 
the transmitted and the incident fields by 

bv(^) = — Im[n(cLi)]L, 

0y{S) = — Re[n(w) — l]L. (3) 

c 

The transmitted field Et results from the interference 
between the incident field Eq and the field scattered in 
the forward direction Es, 


Et — Eq + Eg. (4) 

For effective two-level atoms, we can drop the vecto¬ 
rial nature of the electric fields and represent them as 
scalar quantities. Because of the noninstantaneous re¬ 
sponse time of the medium, the coherent scattered field 
in the forward direction is a continuous function across 
the abrupt change of the incident field. In a FID experi¬ 
ment where the incident field is abruptly switched off at 
t = 0, the intensity of the transmitted field at t = O'*' 
is a direct measurement of the forward scattered inten¬ 
sity in the stationary regime. Its properties are studied 
in detail in mnTj. In particular, the intensity of the 
forward scattering is bounded by 4 times the incident 
intensity (“superflash effect”) HZj. The temporal evo¬ 
lution of the transmitted field, after the abrupt switch 
off of the incident field, is not a simple function having 
only one characteristic decay rate |26j . However, we get 
a clear physical insight by considering only the initial de¬ 
cay time (at t = O'*"), which takes a simple analytical 
expression (see Supplemental Material [35]): 

_ _ /((< = 0+) -/((t = oo) 

dlt/dt{t = 0+) 

2 1-I-exp(—6) — 2exp(—6/2) cos(0) 

r6o(0) 1 — exp(—6/2) cos(6) 


Et{u}) = Eq{(jj) exp 


j)u}L 


where 6 = by{5) and 9 = 9y(S). In Eq. ([^, 6o(0) is the 
optical thickness at resonance and zero velocity. It is 
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linked to 6o(0) by by{0) = bo(0)g(kv/T) where g(x) = 
vAf/S exp (l/ 8 a; 2 ) erfc {l/VSx) /x [26] • 

For small optical thickness, Eq. ([^ reduces to r 5 ( 0 ) = 
g{kv/T)/T at resonance. It is shorter than ro(0) = l/F 
due to the dephasing effect from the motion of the 
atoms. This has already been observed experimentally 
[see Fig. 3(b) of [55] where the transition is Doppler 
broadened, and Fig. 5 of [M] where Doppler broaden¬ 
ing can be ignored]. In our experiments, gikv/V) « 0.16; 
thus, for by(0) = 19(2), we get 6o(0) = 120(10). A direct 
measurement gives a slightly smaller value, &o(0) = 95(5) 
(see Supplemental Material [3S]). The expression of Ty, 
given by Eq. ([^, simplifies to Tc( 0) = 2[6o(0)r]“^ at res¬ 
onance (5 1, 0 = 0) and to r 5 (±oo) = 4[5o(0)r]“^ 

far from resonance (b = 0 and 0 = 0). The solid blue 
curve in Fig. [^is a plot of Ty(S) for kv/F = 3.4. Ty has 
a weak dependence on 6 and v; it depends mainly on 
6 o( 0 ), which can be much larger than the optical thick¬ 
ness 5o(0) seen by a resonant probe at nonzero temper¬ 
ature. This strongly reduces the lifetime of the forward 
scattered field with respect to the atomic lifetime, F”^. 
Equation has a rather simple physical interpretation: 
the second term represents the geometrical properties of 
the propagation inside the medium (change in amplitude 
and phase shift) while the term 2/r5o(0) represents the 
collective behavior of all excited dipoles. It does not de¬ 
pend on the atomic velocity, but only on the atomic den¬ 
sity integrated along the laser direction, because there is 
no Doppler effect for photons scattered in the forward 
direction. Similarly, it does not depend on the detuning 
because all dipoles decay with the same rate F indepen¬ 
dently of the detuning. 

The FID experiment is performed using an AOM as a 
light switching device [see Fig[^a)]. The experimental 
data points, represented by blue open circles in Fig. [^ 
are in reasonable agreement with the theoretical predic¬ 
tion. The evaluation of dlt/dt{t = O'*") is performed on a 
short temporal window (~ 200 ns) after switching off the 
incident probe. While the flash signal has a good signal 
to noise ratio [see Fig. Hb)], the resulting dlt/dt{t — 0 +) 
values from this analysis are noisier. This leads to the 
large statistical errors for Ty. The slight positive sys¬ 
tematic error, also associated to the determination of 
dlt/dt{t = O'*"), comes from the finite response time of our 
experimental scheme, of the order of 40 ns « (500F)“^. 
To check the latter statement, we use Eqs. 0 and 0 to 
numerically compute Eq{ui) in Eq. Q is determined 
from the measured time evolution of the incident inten¬ 
sity. We then apply, on the numerical signal, the same 
procedure used experimentally to extract Ty, resulting in 
an excellent agreement with the experimental data (see 
Fig.§. 

Instead of a FID experiment, we now consider an 
abrupt jump of the phase of the incident field by tt 
[ see Fig. [^c)], at constant incident intensity. The initial 



Figure 2. (color online). Initial decay time of the coherent 
flash versus the probe detuning at kv/V = 3.4. The zero tem¬ 
perature resonant optical thickness is &o(0) = 120. The blue 
open circles and plain curve are, respectively, the experimen¬ 
tal data points and the theoretical curve for an abrupt switch 
off of the probe. The two horizontal black dashed lines give 
the theoretical predictions at resonance and at large detuning. 
The red squares and dashed curve are the experimental data 
points and theoretical prediction for an abrupt phase jump 
of TT. The blue dotted line and the red dash-dotted line are 
numerical predictions taking into account the finite response 
time of the experimental scheme (see text for more details). 


decay time Ty becomes (see Supplemental Material [35]): 


4 1 — exp(—6/2) COS0 

r6o(0) 2 — exp(—6/2) COS0' 


( 6 ) 


We plot this expression as the red dashed line in Fig. [2] 
If the TT phase jump occurs at t = 0, according to Eq. (|^, 
we have Et{t = 0+) = —Eo{t = 0“) -I- Es{t = 0“). To 
observe the largest possible amplitude of the transient 
field, we choose the probe frequency detuning such that 
the interference between Eo(t = 0 “) and Es{t = 0 “) 
is destructive. This condition is necessarily fulfilled 
when the incident field is at resonance. If by{0) ^ 1, 
\Es{t = 0“)| ~ |Ao|, so we expect a coherent flash with 
a peak intensity, A = dJg. The destructive interference 
condition may also happen at a nonzero detuning if the 
phase rotation experienced by Eg is large enough, for ex¬ 
ample if 65 ( 0 ) 3> 1. In our experiment, this situation 
occurs at | 6 | = II.3F (i.e., superflash regime [57]). In 
this context, \Es(t = 0“)| ~ I.SjAoj; thus, the flash has 
a peak intensity A ~ (I -|- 1. 8)^/0 — 7. 8 / 9 . This value is 
slightly below the maximum value 9/o allowed by energy 
conservation, achievable at larger optical thickness. 

The phase jump is performed using an EOM placed 
on the probe laser path [see Fig. [^a)]. The EOM is 
driven by a high voltage controller and has a slew rate 
~ 2 . 37 r rad The two experimental values (red 

squares), corresponding to (5 = 0 and jjj = 11.3F, are 
shown in Fig.[^ They are systematically higher than the 
theoretical prediction for an abrupt phase shift change 
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because of the response time of the EOM driver. Simi¬ 
larly to the FID experiment, we use the experimentally 
measured EOM driver output to numerically compute 
the It{t) signal. The resulting values of the decay time 
(red dash-dotted line in Fig. agree with the experi¬ 
mental ones. 

We now analyze the cooperative emission when a 
square periodic tt phase jump is applied. We observe 
a pulse train with a repetition time Tr [see an exam¬ 
ple in Fig. Hb)] limited by the relaxation time of the 
system. The cooperative emission in the forward direc¬ 
tion dramatically decreases the repetition time below the 
atomic excited state lifetime. 

Bringing the probe on resonance, we plot in Fig. j^a) 
(red dots and solid curve) the intensity contrast Ic of 
the pulse train. We define Ic = max{/t} — (It) as the 
difference between the maximum intensity max{/(} and 
the mean intensity, (/*) = 1 /Tr It{t)dt. We observe 
an excellent agreement between the experiment and the 
theoretical prediction of Eqs. Q and ([^. At long rep¬ 
etition time, i.e., Tr F“^, the system reaches its 
steady state before every phase jump. Hence, we measure 
Ic — 4Jo — {It) — 4/o. We note that {It) — 0 [see the blue 
open circles and dashed curve in Fig. |^a)] and most of 
the incident power is scattered out by single atom fluores¬ 
cence events. In the ^ Jr < F”^ intermediate regime, 
Ic oscillates and can reach a larger value. Moreover, the 
mean intensity {It) rapidly increases to its maximal value, 
Jq. Here, the incident power is almost perfectly trans¬ 
ferred to the pulse train. This interesting result can be 
understood considering cooperativity in forward scatter¬ 
ing. Indeed, its characteristic relaxation time scales like 
[6o(0)r]“^. Therefore, for 6o(0) 3> 1, coherent processes 
relax much faster than single atom fluorescence events. 
The latter are quenched, leading to the good hgure of 
merit at a repetition time shorter than F”^. In other 
words, the emission from the atoms is governed by coop¬ 
erativity. For Jr < Tc, the repetition rate is faster than 
any time scale of the atomic ensemble. Even though the 
probe power is fully transmitted, the contrast Ic tends to 
zero. 

At detuning |(5| = 11.SF, for long repetition time 
(Tr F“^), the pulses have a higher contrast, Ic — 
(1 -I- 1.8)2/o - {It) ^ 7-lJo [see Fig. (b)]. The large 
value of the mean intensity, namely, {It) — O.TIq, is due 
to the small optical thickness, by(S) = 0.4. Hence, most 
of the transmitted power is in a continuous transmission 
mode and not in the pulse train. At the intermediate 
repetition time (t® < Tr < F“^), the pulse contrast and 
the figure of merit are not as good as in the resonant 
case. 

To conclude, we generate pulse trains of short repeti¬ 
tion time using cooperative forward emission in an op¬ 
tically thick scattering medium. We can almost com¬ 
pletely transfer the incident power into the high intensity 
contrast pulse train, quenching the single atom fluores- 
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Figure 3. (color online). Figures of merit of the generated 
pulse train (a) at resonance and (b) at |5| = 11.SF. The red 
solid and blue dashed curves are the theoretical predictions 
for the intensity contrast Ic/Io and transfer efficiency {It)/I q, 
respectively. The red dots and blue open circles are the cor¬ 
responding experimentally measured values. 


cence. This means that, in free space, the cooperativity 
effect can dominate emission from a dilute atomic gas. 
The decay time of the pulses also weakly depends on the 
temperature of the gas and on the probe detuning. An 
interesting extension of this study could be to look for 
quantum signatures in the cooperative emission. 

Finally, we employ the narrow intercombination line of 
strontium as a proof of principle, where the time scales 
are of the order of microseconds. For future practical 
applications, such as a high contrast pulse generator, 
shorter repetition times in the picosecond or subpicosec¬ 
ond regime should be attainable. For this purpose, one 
has to use scattering media with higher optical thickness 
and/or shorter transition lifetime. The fact that cooper¬ 
ativity is robust over thermal dephasing means that we 
can also use a hot vapor of rubidium [6^(0) « 600 at 
110 °C [36]]. We can also use condensed matter systems, 
e.g., a samarium doped fiber [6(0) « 100 |37|], which al¬ 
lows us to bring this technique into the 1.55 fj^im telecom¬ 
munication band. 
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SI 


SUPPLEMENTAL MATERIAL 


Optical thickness measurement 

We employ three different methods to measure the op¬ 
tical thickness. First, we compute the theoretical trans¬ 
mission spectrum for various by{S) values and use these 
profiles to fit the experimentally obtained transmission 
data. This leads to an optical thickness 6^(0) = 19. Sec¬ 
ond, we perform a shadow imaging experiment on the 
—>■ ^Pi broad transition (At, = 461 nm, linewidth 
Ff, = 27r X 32 MHz), where Doppler broadening is negli¬ 
gible. A collimated probe beam with a waist larger than 
the atomic cloud is sent onto the cloud, and the trans¬ 
mission signal It/Io is measured using an electron multi¬ 
plying CCD camera {Andor iXon Ultra 897). Typically, 
the probe frequency is set at a detuning. Si, = 53 MHz, 
to reduce the systematic error in the transmission mea¬ 
surement due to large optical thickness. The optical 
thickness B is computed from the transmission signal, 
B = — log(/t//o), and is related to 6o(0) of the intercom¬ 
bination line by 6o(0) = B (l + dd^/F^) A^/A^. In our ex¬ 
periment, we measure a peak value of 6o(0) = 95(5) using 
this method and a corresponding value of 6 b( 0) = 15(1) 
using by{0) = bo{0)g{kv/T). Third, we carry out shadow 
imaging experiment directly on the intercombination line 
transition. We vary the detuning in a range of 100 kHz 
around the resonance. The value of 5g(0) is deduced us¬ 
ing 


Et{u;) = Eo{u;)U^^, (SI) 


and 

a{ui) 


37rrc3 1 

J-oo '^S-kv + ir/2’ 


(S2) 


which are Eqs. (1) and (2) in the main text. We have 
by{0) = 19(2), a value slightly larger than the one ob¬ 
tained by the second method. 


Initial decay time Ty 


We take t = 0 as the time when the abrupt change 
occurs for the incident field Eq. To calculate the initial 
decay time of the cooperative forward transmitted field, 
we first note that we can rewrite Eq. (5) in the main text 
as 


.(<5) = 


l-\Et{t = oo)\y\Et{t = 0+)\^^ 


rdBt/dt(t=ou\ 

Et{t= 0 +) / 


(S3) 


where Et{t = 0+) = Eo{t = 0+) -I- Es{t = 0“), and 
Et{t = oo) is the steady state transmitted field after the 
abrupt change in the incident field. For the case of abrupt 


extinction, Et{t = 0+) = Es{t = 0“) and Et{t = oo) = 0. 
For abrupt ignition, Et(< = 0+) = Eq and \Et{t = oo)| = 
lEole”^/^. Here, b = by{6) and 6 = 9y{S), the same as 
defined in Eq. (3) of the main text: 

by{6) = —Im[n(a;)]L, 0g((5) = — Re[n(a;) — 1]L. (S4) 


For abrupt phase jump by ip, we have Et{t = O'*') = 
EQe''‘^ + Eg{t = 0“), ignoring the small propagation time 
Ljc in the medium, and \Et{t = oo)| = |Eo|e“^/^. The 
forward scattered field during the steady state regime, 
in both cases of abrupt extinction and phase jump, is 
Eg{t = Q-) = EQe-^/^+^^ -Eq. 


In the denominator of Eq. (S31, we need to com¬ 
pute the time derivative of the transmitted field at t = 
0"*". It can be computed by considering the derivative 
d [i?s(t)e*“*] /dt. The forward scattered field in the time 
domain, Es{t), is related to the incident field in the fre¬ 
quency domain, Eq{oj'), by the following well-behaved 
integral: 


Es{t) = 


^-iUt 



Eo(a;')dw'. 


(S5) 


The integration ranges of the integrals in this Supple¬ 
mental Material, when not specified, are from —oo to oo. 
Eq{u}') is given for the cases of abrupt ignition, abrupt 
extinction and abrupt phase jump of p by: 

i:;o(w') = ^PV^^ + ^^(c^'-a;). (S6) 

27r Lo' — uj 2 


where lo is the frequency of the probe. The Eourier vari¬ 
able corresponding to t is denoted as w'. ^ and rj are —1 
and 1 respectively for abrupt extinction of the probe, and 
e**^ — 1 and 1 -I- e®”^ respectively for abrupt phase jump 
of the probe field. In the case of abrupt probe ignition, 
both ^ and r] are equal to 1. We substitute Eq. (S6) 
in Eq. (S5), noting that the integral involving the Dirac 
delta function goes to zero, to obtain 


_ ^Eq \ - 1 I" /, 

2^ J [ 2c ) 


(S7) 


We work in the regime where S, F, kv loq. For p = 1, 


















S2 


the integral can be evaluated to be 


vanish at t = 0^. Finally, we have 


2tt J 2c 
^Eo &o( 0 ) r 1 


27ri 2 2 y/^% 


dv duj' 




y/^V 


uj' — ujo — kv + iT/2 




for t > 0. In Eq. (S8), 6o(0) = 6tt pc^L/ ujq. The p > 1 
terms, in general, are difficult to evaluate for the general 
time dependence. Nevertheless, at t = O'*', they vanish. 
We take the example of the term p = 2, where essentially 
we have to deal with the following triple integral. 


^ [Esit)e^-^] {t = 0+) = (Sll) 

We then use the fact that Es{t = 0+) = Et{t = 0+) — 
i?o(t = 0+) to obtain 

^(t = 0+) = -eEo^^^-zwEi(t = 0+). (S12) 

dt 4 

Using the above expression, we deduce the initial decay 
time for the case of abrupt probe extinction, 

2 1+ exp(—6) — 2 exp(—6/2) cos(6*) 

r6o(0) 1 — exp(—6/2) cos(6*) 

(S13) 

which is Eq. (5) of the main text. For the case of abrupt 
phase change, we find the initial decay time to be 


dw'dv dU 

(w' — ujo — kv + iT/2){ui' — Wo — kv' + iF/2) 

(S9) 

We rewrite for v ^ v', 


TviS) = 


4 1 — exp(—6/2) COS0 


(S14) 


F6o(0) 2 — exp(—6/2) cos 9' 

This is Eq. (6) of the main text. The initial decay time 
of the flash in the case of abrupt ignition is found to be 


{ui' — ujQ — kv -\- iT/2){ijj' — Wo — kv' + iF/2) 
k{v — v') w' — Wo — fcv + /r 72 

k{v — v') w' — Wo — kv' + zF/2' 

The integration over w' of the above expression can be 
carried out easily, which results in 0 at t = 0+. When 
V = v', we have an integral over a multiple pole of order 
2, which also goes to 0 at t = 0+. Therefore, the term 
with p = 2 is zero at t = O’*". Similar argument can be 
extended to all orders p > 1, showing that all p > 1 terms 


Tv{5) 



(S15) 


We observe again the appearance of the factor 2/F6o(0) 
which arises from the cooperativity among the atomic 
dipoles. 

In the case of an abrupt phase jump, we can further 
choose in the experiment, for p to be equal to the phase of 
Es{t = 0“) relative to Eo{t = 0“). This choice ensures 
a constructive interference after the phase jump. The 
decay time can be simplified to 


To(6) = 


4 |F;,(t = 0+)| 

F6o(0)Uo + |E,(t = 0+)r 


(S16) 






















